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ABSTRACT 

Small  scale  plastic  yielding  at  crack  tips  is  studied  by 
means  of  nonlocal  elasticity.  Plastic  lines  along  the  crack 
line  of  Mode  III  crack  are  modelled  by  an  array  of  disloca¬ 
tions.  It  is  shown  that  plastic  yield  begins  after  a  defi¬ 
nite  value  of  load,  as  a  consequence  of  nonlocality.  The 
length  of  plastic  zone  and  the  dislocation  distrubution  are 
determined  as  functions  of  the  applied  load.  Results  are  in 
good  agreement  with  experimental  observations. 


1 . INTRODUCTION 

One  of  the  fundamental  problems  of  fracture  mechanics  is 
the  determination  of  elasto-plastic  stress  field  near  crack 
tips  in  elastic  solids.  Importance  of  this  knowledge  to 
designers  needs  no  comments.  Traditionally  this  problem  is 
approached  by  means  of  classical  plasticity  theory  which 
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makes  no  reference  to  the  internal  structure  of  materials. 
However,  it  is  well-known  that  the  ultimate  understanding  of 
the  crack  initiation  ,  plastic  yielding  and  fracture  process 
are  intimately  tied  to  the  size  and  distibution  of  defects 
(such  as  impurities,  voids,  dislocations  etc.)  which  exist 
naturally  in  materials.  Consequently,  the  connection  and 
interrelation  of  the  crack-dislocation  assemblies  to  the 
macroscopic  plasticity  and  yielding  are  fundamental  to  the 
fracture  mechanism.  It  is  well-known  that  classical  elas- 
ticty  predicts  an  infinite  stress  at  the  tips  of  a  line 
crack.  This  result  prohibits  the  use  of  the  concept  of 
finite  yield  stress  at  which  plastic  yielding  begins.  As  a 
result  various  ersatz  (such  as  energy,  J-integral,  stress 
intensity  factor,  fracture  toughness  etc.  [1,2]  )  have  been 
invented  to  circumvent  the  difficulty  posed  by  the  stress 
singularity.  Clearly  occurence  of  the  singularity  is  a  def¬ 
inite  sign  of  the  failure  of  classical  elasticity  in  the 
vicinity  of  sharp  crack  tips. 

Fracture  of  a  ductile  solid  with  crack  is  always  accompa¬ 
nied  by  significant  plastic  deformations  in  the  vicinity  of 
crack.  The  plastic  yielding  near  the  crack  tip,  in  the  con¬ 
text  of  classical  continuum  mechanics  was  studied  by  many 
authors  [3-9].  Celebrated  among  these  works  is  the  work  of 
Bilby,  Cottrell  and  Swinden  [5]  who  modelled  the  plastic 
zone  by  an  array  of  dislocations  coplanar  with  the  crack. 


This  model  is  adopted  here  with  basic  departure  being  the 
use  of  nonlocal  continuum  mechanics. 

In  material  science,  still  other  models  are  considered 
relying  on  the  atomic  and  discrete  granular  nature  of 
materials.  Real  materials,  as  distinguished  from  perfect 
crystals,  possess  extremely  complicated  inner  structure. 
Mathematical  analysis  of  crack  tip  problem,  in  the  context 
of  atomic  theory,  faces  extreme  dif f iculities ,  since  the 
complicated,  often  unknown  geometry  of  atomic  distributions, 
large  number  of  dislocations  and  impurities  cannot  be  repre¬ 
sented  with  any  acuracy.  Even  if  this  were  possible  the 
atomic  computations  requires  tedious,  lengthly  and  high  cost 
computations. 

Recently  developed  nonlocal  elasticity  (cf.  10-14)  is 
equipped  with  inner  characteristic  length  mechanism  which 
covers  on  the  one  hand  the  atomic  and  molecular  limits  on 
the  other  classical  elasticity.  Eringen  and  his  co-workers 
have  shown  that  Griffith's  crack  problems  [15-18]  and  point 
dislocations  [19-21]  in  elastic  solids  do  not  possess  singu¬ 
lar  stress  when  they  are  treated  within  the  frame  of  nonlo¬ 
cal  elasticity.  Moreover,  fracture  criteria  based  on  the 
maximum  stress  hypothesis  is  valid  and  it  predicts  correct 
cohesive  stress  for  perfect  crystals.  Several  other  solu¬ 
tions  based  on  the  nonlocal  elasticity  [22-24]  have  shown 
clearly  the  power  and  the  potential  of  the  nonlocal  contin- 
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uum  mechanics. 
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The  present  work  is  concerned  with  the  investigation  of 
the  plastic  yield  at  crack  tip  by  means  nonlocal  elasticty. 
We  consider  a  mode  III  crack  which  containes  a  distributed 
dislocations  of  arbitrary  magnitude  and  lengths  along  the 
crack  line.  Thus  the  plastic  yielding  envisioned  is  small- 
scale.  An  integral  equation  is  obtained  for  the  stress  field 
along  the  crack  line  through  the  condition  that  the  maximum 
stress  reach  the  cohesive  yield  stress.  Computer  solution  of 
this  equation  gives  the  dislocation  density  and  the  plastic 
zone  size  as  a  function  of  applied  load. 

In  section  2  we  present  a  summary  of  basic  equations  of 
linear,  nonlocal  elasticity  for  homogeneous,  isotropic  elas¬ 
tic  solids.  Section  3  contains  Eringen's  result  (25)  on  the 
stress  field  due  to  interaction  of  single  dislocation  with 
crack.  These  results  constitute  the  Green's  function  in  the 
formulation  of  the  distributed  dislocations.  Section  4  is 
devoted  to  the  solution  of  the  problem  and  the  results 
obtained  by  computer  work.  Last  section  gives  a  discussion 
and  physical  significeuice  of  these  calculations. 

2. BASIC  EQUATIONS 

In  several  previous  paper  (for  exaunple  cf.  [26]),  Eringen 
has  shown  that  the  field  equations  of  linear,  nonlocal  elas¬ 
ticity  for  homogeneous,  isotropic  solid,  with  vanishing  body 
forces,  consist  of  the  equations  of  equil ibrium 


t. .  • (X)  =  0 

ij, 


(1) 


and  the  stress-strain  relations 


^ij  ^2)=b‘^«(  Ix-x'  I  (2^)civ(x'  ) 

where  is  the  classical  Hookean  stress 

0ij(x)=Xe^^(x)Cl.t2vei.(x) 

which  is  expressed  as  a  linear  function  of  the  strain 
2ei.(x)  =  u.^j(x)  ^  u.  .U)  . 


(2) 


(3) 


(4) 


Here  u^(x)  is  the  displacement  vector,  X  and  v  are  the  Lame' 
elastic  constants  and  indices  following  comma  denote  the 
partial  differentiation  with  respect  to  space  variable. 
®  (  li5”2S' I  /  E )  is  the  "attenuation  function"  which  brings  the 
effect  of  strains  at  distant  points  in  body  to  the  stress 
field  at  a  reference  point  x  .  This  function  depends  on  a 
characteristic  length  e  .  The  integral  in  (2)  is  over  volume 
of  the  body  denoted  by  B.  The  function  a  is  a  decreasing 
function  of  I2-2' |  .  Since  intermolecular  forces  die  out 
fast  with  distance,  the  effective  range  of  o  is  of  the  size 
of  molecular  order.  For  perfect  crystals  it  is  of  order  of 
the  lattice  parameter  while  for  granular  and  porous  solids 


it  is  of  the  order  of  micron.  If  we  denote  the  external 
characteristic  length  by  L  and  an  internal  characteristic 
length  at  natural  state  of  the  body  by  t  then  the  ratio 

e  =  eot/L  =  e/L  (5) 

where  e,  is  a  dimensionless  material  constant,  determines 
the  range  of  validity  of  the  continuiim  theory.  For  e=0  we 
have  the  classical  (local)  elasticity  and  for  e?to  nonlocal 
elasticity  capable  in  dealing  with  microscopic  and  atomic 
phenomena.  In  fact  it  can  be  shown  that  for  an  appropriate 
function  o  the  nonlocal  elasticity  gives  exactly  the  same 
stress  field  at  atomic  points  (cf.  Eringen  [27]).  This  func¬ 
tion  can  be  determined  in  various  ways:  by  experiment,  by 
statistical  averages  applied  to  atomic  theories  and  by  com¬ 
paring  dispersion  relations  of  elastic  waves  with  the  phonon 
dispersion  curves.  For  example,  Ari  and  Eringen  [18]  using  a 
two  dimensional  lattice,  have  obtained  an  excellent  match 
in  the  entire  Brillouin  zone,  with  phonon  dispersion  curves, 
for 

o(  Ixl  )  =  (2ire*)-^K^(»^T^/E)  (6) 

where  is  the  modified  Bessel’s  function  of  first  kind, 
o 

For  other  kernels  in  1,  2  and  3  dimensions  see  [24]. 


It  is  interesting  to  note  that  the  expression  given  by 


(6)  is  the  Green's  function  of  the  linear  operator  L=1-e*V* 
i  .  e .  , 

(l-E*7Ma(l2-2'l)=6(|x-£'l)  (7) 

where  6(x)  is  the  Dirac  delta  measure.  This  feature  of  the 
function  allows  us  to  invert  the  constitutive  equation 
(2)  of  nonlocal  elasticity: 

(l-E^V*)t.  .=0.  .  (8) 

With  this  apparatus  at  hand,  Eringen  [25]  gave  the  solution 
of  the  boundary  value  problem  of  nonlocal  elasticity  on  the 
interaction  of  a  dislocation  with  crack,  relevant  to  the 
present  work.  This  is  summarized  in  the  following  section. 

3. FORMULATION  OF  THE  PROBLEM 

The  main  purpose  of  this  study  is  to  describe  the  plastic 
zones  at  the  tips  of  a  crack  for  Mod  III  problem.  The  plas¬ 
tic  zones  will  be  modelled  by  an  array  of  dislocations.  We 
consider  a  homogeneous,  isotropic,  elastic  medium  of  infi¬ 
nite  extend  weekend  by  a  crack  located  at  -c  S  x^  S  c  , 
X2=0  ,  —  <  x^  <  -  and  is  subject  to  a  constant  antiplane 

shear  load  ip  at  X2=+*  We  assume  that  there  exist  a  screw 
dislocation  whose  Burger's  vector  b  is  parallel  to  the  x^ 


ly;:.' 


axis  and  intersects  the  plane  *^=0  at  the  point  S(  , 

X2=0).,(see  Fig.l)  .  The  classical  solution  of  this  prob¬ 


lem  (i.e.  interaction  of  a  crack  and  a  dislocation)  has  been 


given  by  Louat[5]  in  a  general  form.  The  classical  solution 


of  this  problem  has  also  been  given  by  Eringen[25]  in  a  more 


appropriate  form  which  may  be  expressed  as  follows: 


023(x^,X2)-io.,.^(XT  ,x,)  = 


(z*-c"r'^{ToZ+(2itrVbIl+(5*-c*)'^(2-?)”M> 


where  i=v'-l  ,  z=x^  +  iX2  is  the  complex  variable  ,  z=x^-iX2 


being  its  complex  conjugate.  The  net  content  of  disloca¬ 


tions  in  the  crack  is  assumed  to  be  zero.  This  expression 


contains  four  different  response  of  the  medium: 


(i)  Stress  field  due  to  applied  load  <>23 


a23(Xi,X2)-iOi3(Xi,X2)=t, 


(ii)  Stress  field  due  to  the  response  of  the  crack  to  the 


applied  load: 


O23  ( X ,  X2 ) -io ( Xj^ ,  X2  )=T  0 { z  ( z* -c*  )‘V’* - 1 } 


(iii)  Stress  field  due  to  the  existence  of  the  disloca¬ 


tion  at  x^=5  ,  X2=0: 


023(Xi,X2)-iOi3(Xi,X2)  =  (2Tr)-*vb(z-5)‘‘ 


(iv)  Stress  field  due  to  the  interaction  of  the  disloca¬ 
tion  with  the  crack: 

( 2it r yb{  ( Z"-c*  [  1+  (2-$ )-‘  ]  -  (2-^  T >  ( 13  ) 

It  is  clear  that  the  expression  (9)  has  square  root  singu¬ 
larities  at  the  crack  tips  and  1/r  singularity  at  the  disla- 
cation.  For  the  nonlocal  solution  of  this  problem  Eringen 
[25]  has  argued  that  the  nonlocal  stress  field  does  not 
possess  any  singularity.  To  obtain  the  nonlocal  stress  field 
it  was  necessary  to  find  the  solutions  of  the  following  par¬ 
tial  differential  equation 

(l-c*7Mt  =  0  (14) 


where 


t(Xi,X2)  t22(x^,X2)-it^2(^l'*2^ 

and  0  is  given  by  (9).  It  is  easy  to  verify  that  t=o  is  a 
particular  solution  of  Eq.(14).  The  general  solution  of  the 
homogeneous  equation  [(14)  with  o=0]  that  possesses  proper 
symmetry  regulations  and  vanishes  at  infinity  can  be  formed 
from  the  generic  solution 

t  =  K^(r/E){A^e^''®+B^e'^''®) 


(16) 


where  A  ,  B  and  v  are  constants  ,  K  (p)  is  the  modified 
Bessel's  function  of  first  kind  and  (r,0)  are  the  plane 
polar  coordinates. Considering  the  limiting  behaviour  of  mod¬ 
ified  Bessel's  functions  of  first  kind  with  integer  order, 
for  small  values  of  argument  we  have  (c.f.  (9.6.9)  in  [28]) 

(17) 

K^(z)  ~  .5r(v)(  .Sz)"'' 

and  the  special  representation  of  (c.f.  (10.2.17)  in 

[28]) 

(18) 


The  nonlocal  stress  field  is  then  given  by  [25] 


where 


and 


_ V 


-riA  ^-i0i/2 


t23(x^,X2)-itj^2(’^l'^2^“  C^(TrE/2r^r^e"  e 


C2K^(r3/c)e^^®’"®‘^  +  ( r^r2  T'^e^  ^  ^  {tore"^® 


+  (2ir)-‘yb{l+[5(5+2c)]V*  (  r^e^® ‘-5  )-^ } 


-Vi 


C^=-(c/irE)  {T,  +  (2Trc)-‘vb[l-(l+2cA)^]> 


C2=-  (2iie  r‘yb 


(19) 


(20) 


(21) 
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For  the  purposes  of  this  study  we  need  to  consider  an  anti¬ 
symmetric  dislocation  distribution.  Consequently  we  super¬ 
pose  the  stress  field  due  to  the  dislocation  located  at 


(X^=-C  ,  ^2“°^  with  -b  Burger's  vector  to  (19).  With 
this  the  singularity  is  removed  from  the  crack  tip  at  x=-c  . 
The  stress  distribution  may  be  expressed  as 

(c/2r^)^{to-(ircrwb[U5+2c)  rV^(C+c)}e‘^i/^ 


(c/2r2)V^{to  +  (Trc)-‘yb[U^+2c)  (C+c ) 
(2irE)-Vb{e^®^K^(r3/E)+e^®'*K^(r4A)}  (22) 


where  9^  and  r^  are  the  polar  coordinates  as  shown  in  Fig. 2 
.  Since  we  are  dealing  with  small  scale  yielding  near  crack 


tip  we  need  only  the  t23  component  of  stress  on  the  plane 
coplanar  with  the  crack.  From  (22)  we  have 

t23(x,0)=t.^(x)+t'''^(x) 
t^(x)=-[o{  [x(x+2c)  fV^(x+c)-(c/2x) 
t'='^(x,C)  =  (if^W) 

{[x(x+2c)  rV*[C(C+2c)  ]  v4(x-U‘‘(x+^+2c)-‘(?+c)  + 

(2c)-V^[U5+2c)]*V^(5+c)x-V*e"*/''- 
(2Er[sgn(x-UK3(  |x-S|/e  )+K^(  (x+C+2c)/e  )  ] } 

(23) 


where  x  denotes  the  distance  from  the  crack  tip  at  right 
hand  side.  The  stress  field,  due  to  an  array  of  disloca¬ 
tions  in  the  interval  (a,b)  with  density  b(C)  can  now  be 
expressed  as: 

t(x)  =  t^(x)  +  t‘='^(x,Ub(?)dC  (24) 


4. SOLUTION  OF  THE  PROBLEM 

From  the  point  of  our  intrest  we  aim  to  solve  the  Eg. (24) 
for  t(x)=ty,  d<x<e ,  where  t^  is  the  cohesive  stress  between 
atoms  and  d,e  are  the  coordinates  of  the  plastic  zone.  The 
equation  which  we  need  to  solve  is  a  Fredholm  integral  equa¬ 
tion  of  first  kind  with  a  smooth  kernel.  The  well-known 
ill-posed  character  of  such  equations  increases  with  the 
smoothness  of  the  kernel.  In  general,  all  information  pro¬ 
vided  for  these  equations,  either  analytical  or  numerical 
point  of  view,  is  crude  rule  of  thumb  rather  than  a  precise 
guide  to  the  nature  of  the  problem.  For  analytical  and 
numerical  treatment  of  Fredholm  integral  equation  of  first 
kind  [29-32]  can  be  referred. 

Among  the  numerical  techniques  for  solving  the  Fredholm 
integral  equation  of  first  kind,  the  so-called  Minimal  Least 
Square  Solution  is  perhaps  the  most  reliable  and  handy  one. 
According  to  this  method,  the  unknown  function  is  approxi¬ 
mated  by  a  set  of  basis  function 


b(n  =  Z  b,  0,(5) 


and  the  unknown  coefficients  b.  are  sc  described  that  the 

1 

difference  between  the  results  obtained  from  the  integral 
with  (25)  and  the  desired  result  will  be  minimum  with 
respect  to  a  norm 


I 

I 


n 

r  b, 
i  =  l  ^ 


(Ud^  +  t‘^(x)-ty|  j-^min 


d<x<e 


(26) 


In  some  cases,  some  restrictions  on  b^  may  accompany  to  this 
minimization  problem. 

As  a  first  approach 

^.(5)  =1  q  S  ^  (27) 

can  be  chosen.  Here  C-  and  denote  the  coordinates  of  a 

subinterval  in  (a,b).  With  this  set  of  basis  function  it  can 

cd, 

be  constructed  a  simple  algorithm  since  t  (x,0  can  be 
integrated  with  respect  to  i  : 


ti(x)  =  I  t''^(x,UdC 


(28) 


Before  performing  the  integration  let  us  normalize  the  dis¬ 
tances  with  the  half  crack  length  c  and  the  stresses  with 
the  cohesive  stress  t 


z  =  x/c 


T1  =  K/c 


R  =  T  /t 
o  y 


(29) 


In  terms  of  these  new  variable  we  have 


t‘^(z)=R{  [z(z+2)  r^(z  +  l)-(2zrV2e“^‘^/'^^^} 
t^^(Z,Tl)  =  (TTC)‘^y 

{[z(z+2)  rV*[T,(T)+2)  ]  Vi*(z-Ti)-Mz+Ti+2)->(n  +  l)  + 

(l//z)[Ti(n  +  2)r^(Tl  +  l)z-V^e“^^/^^^- 

(c/2e  ){sgn(Z“ii  )K^  (  (c/e  )  1  z-n  |  1+K^  [  (c/e  )  (x+ti+2  )  ] }} 

(30) 

To  calculate  the  integral  (28)  we  consider  the  following 
indefinite  integrals: 

/[x(x+2)  ]  v!*  (y-x)-^dx=-[x(x+2)  ]  '^+[y(y+2)  ]  ^ 

{ln{  [y(y+2)  ]  V<*(x(x+2)  ]  ^-(y+l)(y-x)+y(y+2)>- 
ln(y-x)}+(z+l)ln{(x(x+2) ] + ( y-x)- (y+1 ) }  (31) 


It  should  be  noted  that  the  definite  value  of  this  integral 
exists  in  the  sense  of  Cauchy  Principle  Value  if  z  takes 
place  in  the  integration  interval. 


;[x(x+2)  1  Vi^(y+x+2)'^dx=-lx(x+2)  ]*A+[y(y+2)  ] ’A 
{ln{  [y(y+2)  ]  ‘^[x(x+2)  ]  ^-(y+l)(y+x+2)+y(y+2)}- 
ln(y+x+2)}+(z+l)ln{  [x(x+2) l  +  (y+x+2)-(y+l))  (32) 


;  [x(x+2)  rV^(x+l)dx=(x(x+2)  )  ^ 


(33) 


/sgn(y-x)KT ( ly-xi )dx=K  ( ly-x| ) 


(34) 


Also  the  definite  value  of  this  integral  exists  in  the  sense 
of  Cauchy  Principle  Value  if  z  is  in  the  integration  inter¬ 
val.  With  the  aid  of  these  results  can  be  expressed 

as  follows: 

(z)  =  (^-'y){(2/A)(R2-R^)  +  lnE^  +  lnE2  + 


(l/2)[K^((c/E)!A^i )-K^((c/E)iA2l )- 

(  (  c /E  )  ( B+B^  )  )  (  (  c/E  )  ( B+B2))] 


where 


7i^=c-Ma+(b-a)  (i-l)/N]  i=l,2,  .  .  .  ,N+1 


A=[z(z+2)]'^  ,  +  l 

B=z  +  1  ,  +  '  ®2~'^i'^^ 

E^=|A2/A^|  (B+B^rCB+B^) 

.2, 


E-t=(A.Rt-B.At+A  ) (A.R,-B.B,-1) 


E22=(A.R2-B.A2+A^) (A.R2-B.B2-I) 


^2  ^21^^22 


Let  us  define  the  residuals 


n 


:.=  Zb.  'i'.(z.)  +  t  (z.)-t 
1  1  1'  '3'  j 


d<Zj  <e 


(35) 


(36) 


(37) 


Now,  we  wish  to  describe  the  unknown  coefficients  b^  by  min¬ 
imizing  the  squares  of  residuals 


m 


Vi’4  4’4V4  ••i  i'l  l*!  * 


F 


(38) 


m 


=  Z 

j  =  l 


(r.) 


2 


It  should  be  noted  that  a  necessary  condition  for  the  exis¬ 
tence  of  a  set  of  ' s  which  minimize  the  functional  F  is 
that  m  >  n  .  The  fundamentals  of  the  solution  techniques  of 
these  kind  of  problems  Dennis  and  Schnabel  [33]  can  be  con¬ 
sulted.  Since  the  experimental  results  [34]  indicate  that 
the  dislocations  near  a  creack  tip  have  always  the  same  sign 
we  solve  the  Eq.(38)  with  the  constrains 

b^  >  0  (39) 

In  the  solution  of  this  problem  the  following  points  are 
important : 

(  i)  The  begining  and  end  points  of  the  dislocation  zone 
(a,b)  and  of  the  plastic  zone  (d,e)  are  the  unknowns  of  the 
problem  for  a  given  applied  load.  There  is  no  appearent 
mathematical  reason  to  take  a=d  ,  b=e  which  is  the  case  in 
classical  approach.  Neverthless,  good  results  (in  stress 
and  dislocation  distribution)  are  obtained  when  they  are 
chosen  very  close  to  each  other. 

(ii)  An  important  difference  between  classical  and  nonlo¬ 
cal  approach  is  the  occurence  of  the  plastic  yield  at  crack 
tip.  As  is  well  known,  in  classical  approach,  plastic  yield¬ 
ing  occures  for  every  value  of  applied  load.  But  this  is  not 
true  in  nonlocal  approach.  Since  the  response  of  a  crack  to 


the  external  load  possesses  an  extremum  near  the  crack  tip 
in  nonlocal  approach  the  plastic  yielding  will  begin  after  a 
definite  value  of  applied  load  for  which  the  extremum  value 
of  elastic  response  exceeds  a  prescribed  limit  (maximum 
stress  hypothesis).  So,  the  nonlocal  elasticity  enables  us 
to  define  a  critical  value  of  applied  load  after  which  the 
plastic  yielding  occures  at  the  crack  tips. 

(iii)  The  coordinates  of  dislocation  zone  (or  plastic 
zone,  assuming  the  coordinates  of  both  zones  are  Seune)  are 
unknown.  An  examination  of  the  behaviour  of  the  nonlocal 
elastic  stress  field  near  crack  tip  clearly  indicates  that 
the  beginning  point  of  dislocation  zone  is  very  close  to  the 
point  at  which  the  elastic  stress  distribution  reaches  its 
extremum  value.  The  end  point  of  dislocation  zone  is  a  func¬ 
tion  of  the  applied  load.  To  determine  this  point  for  a  giv¬ 
en  value  of  applied  load  it  is  necessary  to  solve  a  diffi¬ 
cult  nonlinear  problem  in  which  applied  load  is  unknown. 
Instead  of  tedious  computations  required  by  this  process  we 
choose  the  applied  load  as  unknown  for  a  given  value  of  the 
end  coordinate  of  plastic  zone.  This  approach  provides  much 
simplicity  and  economy  in  the  volume  of  computations. 


The  problem  to  be  solved  is  to  determine  the  dislocation 
densities  b^  and  the  applied  load  Tq  so  that  the  functional 
F  given  by  (38)  reaches  its  minimum  value  under  the  con¬ 
strains  (39).  The  main  parameters  of  the  problem  are  the 
internal  characteristic  length  (t)  and  the  half  crack  length 
(c).  The  results  obtained  for  e  =  .00001  (cm.)  and  c  = 
.001(cm.)  are  given  in  a  series  of  figures  (Fig. 3-10).  In 
the  figures  marked  with  (a)  are  shown  the  stress  distribu¬ 
tion.  The  curves  marked  as  ’’elastic"  and  "interaction"  show 

the  contribution  of  the  term  t''  and  t^^  which  are  given  by 
(30),  respectively.  The  curves  marked  as  "elasto-  plastic" 

show  the  totality  of  these  terms.  As  is  seen  clearly  from 
these  figures  the  behaviour  of  elasto-plastic  stress  field 
is  quite  reasonable;  stress  increases  to  a  definite  value 
(ty),  stays  unchanged  (plastic  zone)  and  decreases  to 
another  value  (applied  load)  asymptotically.  The  increasing 
and  decreasing  parts  of  these  curves  are  the  elastic 
response  of  the  medium.  Another  important  aspect  is  "dislo¬ 
cation  free  zone"  which  is  clearly  seen  in  the  figures  mark¬ 
ed  with  (b)  in  which  the  dislocation  distribution  are  shown. 
Dislocation  free  zone  which  is  observed  by  in-situ  electron 
microscope  experiments  [34]  is  a  natural  result  of  nonlocal 
approach  i.e.  needs  no  further  assumption.  In  another 
series  of  figures  (Fig. 11-13)  the  results  obtained  for  e  = 


.00001(cm.)  and  c  =  .l(cm.)  are  given.  The  plastic  zone  is 
quite  small  as  compared  with  the  half  crack  length,  in  these 
figures.  For  larger  plastic  zone  it  is  neccesary  to  take 
large  number  of  division  of  the  dislocation  zone  which  makes 
computation  cost  high.  If  the  division  number  is  not  suffi¬ 
ciently  large  then  approximation  is  not  good  enough  as  can 
be  observed  from  Fig.  13  .  In  the  last  figure  (Fig.  14)  the 

end  coordinate  of  plastic  zone  obtained  by  nonlocal  and 
classical  approach  are  shown.  The  nonlocal  result  is  remar- 
kabely  smaller  than  the  classical  result.  But  this  is  quite 
reasonable  because  there  is  no  plastic  yielding  for  some 
values  of  applied  load  smaler  than  a  definite  limit  in  non¬ 
local  approach  while  classical  approach  gives  some  plastic 
zone  length. 

REFERENCES 

[IJ.A.A.  Griffith,  Phil.  Trans.  Roy.  Soc .  London,  Ser.A 
221,163,1921. 

[2] .G.C.  Sih  and  H.  Liebowitz,  Mathematical  Theories  of 
Brittle  Fracture  (in  Fracture  Vol.2,  ed.  H.  Liebowitz),  Aca¬ 
demic  Press, New  York  ,1968. 

[3] .D.S.  Dugdale,  J.  Mech.  Phys.  Solids,  8,  100,  1960. 

[4] .  G.I.  Barrenblatt,  Adv.  Appl.  Mech.,  7,  55-129,  1962. 


[5J.B.A.  Bilby,  A.H.  Cottrell  and  K.H.  Swinden,  Proc. 
Roy.  Soc .  London,  Ser.  A272,  304,  1963. 

[6] .N.  Louat,  The  Interaction  of  Crack  and  Dislocations, 
Proceedings  of  the  First  Intenational  Conference  on  Frac¬ 
ture,  Sendai,  Japan,  117,  1965. 

[7] .C.  Atkinson,  The  Interaction  Between  a  Crack  and  a 
Dislocation,  Int.  J.  Fracture  Mech. ,  2,  567,  1966. 

[8] .J.P.  Hirth  and  J.  Lothe,  Theory  of  Dislocations, 
McGraw  Hill,  New  York,  1968. 

[9] .R.W.  Lardner,  Mathematical  Theory  of  Dislocations  and 
Fracture,  University  of  Toronto  Press,  Toronto,  1974. 

[10] . A. C.  Eringen,  A  Unified  Theory  of  Thermomechanical 
Materials,  Int.  J.  Engng.  Sci.,  4,  179,  1966, 

[llj.E.  Kroner,  Continuum  Mechanics  and  Range  of  Atomic 
Cohesion  Forces,  Proceedings  of  the  First  International  Con¬ 
ference  on  Fracture,  Sendai,  Japan,  27,  1965. 

[12].D.  Rogula(ed.),  Nonlocal  Theory  of  Material  Media, 
Springer  Verlag,  1982. 

[ 13 j . I .A.Kunin,Elastic  Media  with  Microstructure:  Vol.I 
and  II,  Springer  Verlag,  1983. 

[14] .A.C.  Eringen,  D.G.B.  Edelen,  On  Nonlocal  Elasticity 
,  Int.  Journ.  Engng.  Sci.,  10,  233,  1972. 

[15] .A.C.  Eringen,  C.G.  Speziale  and  B.S.  Kim,  Crack  Tip 
Problems  in  Nonlocal  Elasticity,  J.  Mech.  Phys .  Solids,  25, 
339,  1977. 


[16J.A.C.  Eringen  ,  Line  Crack  Subject  to  Shear,  Int. 
Journ.  Fracture,  14,  367,  1978 

[17] .A.C.  Eringen,  Line  Crack  Subject  to  Antiplane  Shear, 
Fracture  Mechanics,  12,  211,  1979 

[18] . N.  Ari  and  A.C.  Eringen,  Nonlocal  Stress  Field  at 
Griffith  Crack,  Crys.  Latt.  Def.  and  Amorph.  Mat.,  10,  33, 
1983. 

[19] . A.C.  Eringen,  Edge  Dislocation  in  Nonlocal  Elastici¬ 
ty,  Int.  J.  Engng.  Sci.,  15,  177,  1977. 

[20] . A.C.  Eringen,  Screw  Dislocation  in  Nonlocal  Elastic¬ 
ity,  J.  Phys.  D:AppL.  Phys.,  10,  671,  1977. 

[21] . A.C.  Eringen  and  F.  Balta,  Screw  Dislocation  in  Non¬ 
local  Hexagonal  Crystals,  Crystal  Lattice  Defects,  7,  183, 

1978. 

[22] . A.C.  Eringen,  Linear  Theory  of  Nonlocal  Elasticity 
and  Dispersion  of  Plane  Waves,  Int.  Journ.  Engng.  Sci.,  10, 
425,  1972. 

[23] . A.C.  Eringen,  On  Rayleigh  Surface  Waves  with  Small 
Wave  Length,  Letters  in  Appl.  Engng.  Sci.,  1,  11,  1973. 

[24] . A.C.  Eringen,  On  differential  Equations  of  Nonlocal 
Elasticty  and  Solutions  of  Screw  Dislocation  and  Surface 
Waves,  J.  Appl.  Phys.,  54(9),  4703,  1983. 

[25] . A.C.  Eringen, Interaction  of  a  Dislocation  with  a 
Crack,  J.  Appl.  Phys.,  54(12),  6811,  1983. 


[26] . A. C.  Eringen, Nonlocal  Elasticity  and  Waves,  Contin¬ 
uum  Mechanics  Aspects  of  Geodynamics  and  Rock  Fracture 
Mechanics,  (ed.  Thoft-Christensen) ,  D.  Reidel  Publishing 
Co.,  Dordrecht-Holland,  81-105,  1974. 

[27] . A. C.  Eringen,  Theory  of  Nonlocal  Elasticity  and  Some 
Applications,  Res.  Mechanica,  21,  313-342,  1987. 

[28] . M.  Abramowitz  and  I. A.  Stegun  (editors).  Handbook  of 
Mathematical  Functions,  National  Bureau  of  Standarts 
Applied  Mathematics  Series. 55  Washington,  1970. 

[29] .C.T.H.  Baker,  The  Numerical  Treatment  of  Integral 
Equations,  Clarendon  Press,  Oxford  ,1977. 

[30] .F.G.  Tricomi,  Integral  Equations,  Interscience  Pub¬ 
lishers,  Inc.,  New  York  1957. 

[31] .C.T.H.  Baker  and  G.F.  Miller,  Treatment  of  Integral 
Equations  by  Numerical  Methods,  Academic  Press,  London 
, 1982 . 

[32] .L.M.  Delves  and  J.L.  Mohamed,  Computational  Methods 
for  Integral  Equations,  Cambridge  University  Press,  1985. 

[33] .J.E.  Dennis  and  R.B.  Schnabel,  Numerical  Methods  for 
Uncostrained  Optimization  and  Nonlinear  Equations,  Prentice- 
Hall,  Inc.,  Englewood  Cliffs,  New  Jersey,  1983. 

[34] .S.M.  Ohr,  J.A.  Horton  and  S.J.  Chang,  Direct  Obser¬ 
vation  of  Crack  Tip  Dislocation  Behaviour  During  Tensile  and 
Cyclic  Deformation,  Proc.  Second  Int.  Sym.  on  Defects,  Frac¬ 
ture  and  Fatigue,  Martinus  Nijhoff  Publishers,  The  Hauge, 


!Sl 


■;S 

'll'!* 

Ll 


I 

O  u 
^  o 


®  2 
Ebt 


(  ^'\/'\  )SS0a'^S  I«U0ISU9UITa-U0K 


Dimensional  Distance 


dimensional  Distance  From  Crack  Tip  (  x/c ) 


dimensional  Distance  From  Crac! 


Non-Dimensional  Distance  From  Crack  Tip  (  x/c ) 


dimensional  Distance  From  Crack  Tip  (  x/c ) 


Non-Dimensional  Distance  From  Crack  Tip  (  x/c ) 


elastic 


Non-dimensional  Distance  From  Crack  Tip  (  x/c ) 


Non-dimensional  Distance  From  Crack  Tip  (  x/c ) 


fwr^  suoi^isooisTQ  10  Jdquiniil 


cation  Distribution  Near  Crac 


Dislocation  Distribution  Near  Crack  Tip 


Dimensional  Distance  From  Crack  Tip  (  x/c ) 


0.001  0.002  0.003  0.004 

Dimensional  Distance  From  Crack  Tip  (  x/'c  1 


Non-dimensional  Distance  From  Crack  Tip  (  x/c 


Non-dimensional  Applied  Load 


